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Abstract 

A statistical model of loops on the three-dimensional lattice is proposed and 
is investigated. It is 0(n)-type but has loop fugacity that depends on global 
three-dimensional shapes of loops in a particular fashion. It is shown that, 
despite this non-locality and the dimensionality, a layer-to-layer transfer matrix 
can be constructed as a product of local vertex weights for infinitely many points 
in the parameter space. Using this transfer matrix, the site entropy is estimated 
numerically in the fully packed limit. 
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1 Introduction 



Loop models are interesting examples of statistical models of extended objects. 
They are related to the O(n) spin model [1, 2], a surface growth model [3], the 
self-avoiding walk [4], the protein folding problem [5], and so on. It includes 
the fully packed loop model [6] and the Hamiltonian cycle problem [7, 8, 9] as 
particular limits. 

The partition function of an O(n) loop model on a lattice with ./V sites at 
the inverse temperature x is given by 

Z ioop (n,x- 1 )=Y / X Ars{c) ~ N n^ (c) (n,x E M). (1) 

cGC 

The summation is taken over the set C of all the non-intersecting loop configu- 
rations drawn along links of the lattice. The number of loops and that of sites 
visited by them are denoted by A/l(c) and A/s(c), respectively. 

One may hope to study the model (1) by the transfer matrix approach. For 
n G Z + , this is done in a simple way; one introduces link variables whose values 
are either occupied states with one of n colors or an unoccupied state and lets 
them interact on sites. A transfer matrix is written as a product of vertex 
weights straightforwardly. 

For n ^ Z_|_, however, the partition sum (1) cannot be rewritten in terms 
of local degrees of freedom such as link variables in a simple way. It is not 
trivial to have a local transfer matrix 1 . I say a transfer matrix is local when its 
component is written as a product of weights each of which is determined by 
the local state configuration around a lattice site. 

It is surprising that, in two dimensions, n $jL Z + models admit a mapping 
onto a state sum model with a local vertex weight and thus have local trans- 
fer matrices [10, 11, 12]. In fact, by choosing s e C satisfying n = s + s _1 , 
Z\oop( n , x" 1 ) can be written as 

^oo P (n,x- 1 )=^x^( c )- JV ( S + S - 1 )^( c )=^x^( c )- JV J] s±\ (2) 
c ec cec Lec(c) 

where C is the set of loop configurations with a direction associated with each 
loop. The set C(c) consists of all the directed loops in a configuration c. A loop 
with the (counter-)clockwise direction is given a weight s +1 (s _1 ). This weight 
can be realized by associating s +1 / 4 (s -1 / 4 ) with each right(left)-turn site and 
the model can be regarded as a state sum model with a local vertex weight. 
This trick has made the study of two-dimensional loop models very fruitful. 

Physics of loops in three dimensions is very attractive. It is realistic in 
the context of condensed matter physics. There has been a continuous suspi- 
cion that two-dimensional ones have missed some important ingredient in real 
physics, e.g. the protein folding problem. Three-dimensional loops have also 
rich mathematical structures. For instance, loops can be knotted or linked in 

lr The use of the connectivity basis is discussed in subsection 5.4. 
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three dimensions [13]. It is noted that a number of attractive proposals have 
been made to generalize the loop model to higher dimensions [14, 15]. 

The analysis of loop models and their generalizations in higher dimensions is, 
however, extremely hard to perform. Needless to say, the number of configura- 
tions increases considerably. For fugacity n Z + , which includes the interesting 
case of the self-avoiding walk (n = 0) , no way of constructing local transfer ma- 
trices is known. This is because specialties of two dimensions cannot be used to 
simplify problems any more. The mapping (2) makes use of the fact that the 
a directed loop in two dimensions turns around just once either clockwisely or 
counter-clockwisely. It appears that this kind of tricks never works in higher 
dimensions. 

In this article, I propose a model which generalizes (1) in a fashion specific to 
three dimensions. It is furnished with loop fugacity that depends on the global 
three-dimensional shape of loops. I show that, despite this generalization which 
makes the model even more non-local, a local transfer matrix for the system 
can be constructed for a number of choices of fugacity. These choices include 
the ones that give zero or non-integer weight to loops. 

This article is organized as follows. In section 2, I define a loop model in 
three dimensions generalizing (1). Its local transfer matrix is constructed for a 
family of points in the parameter space in section 3. In section 4, this transfer 
matrix is numerically diagonalized to yield an estimate of the site entropy in 
the fully packed limit = 0. In the last section 5, I discuss my results and 
their relation to combinatorial problems. In an appendix, a technical issue on 
the block-diagonalization of the transfer matrix is addressed. 

2 Generalized fugacity 

I define a statistical model of loops on the three-dimensional simple cubic lattice 
Z 3 = {Y^i=i TUiei € R 3 \mi e Z} , e, • ej = 5ij. The partition function is given 
by 

^ 1 oop[n](.T- 1 )=^^( c )- JV J] n(A(L)), (3) 
cec Lec(c) 

where C(c) is the set of loops in a configuration c. The loop fugacity n is now 
promoted to a function which depends on the shape of L e £(c) through a 
quantity A(L) G R defined below. 

To define A(L), one begins with associating a closed trajectory on the unit 
sphere with each loop L. One picks a direction for L. On every point xeLcR 3 
except for sites where L makes a turn, there is a unit tangential vector v(x) 
to L; it is either of ±ej, % = 1,2, 3. One may regard v(x) as a mapping from 
L \ ('turn-sites') to the unit sphere S 2 . 

As one walks along L, v(x) jumps from a point to another on S 2 . One can 
naturally interpolate these points to define a continuous trajectory v : L — > S 2 . 
One has only to declare that v(x) moves along the geodesic (of length ^n) on 
S 2 at each turn-site. This is equivalent with smoothing a loop in neighborhoods 
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(a) (b) (c) 



Figure 1: The definition of A(L). (a) The original loop L with a direction 
associated, (b) The loop smoothed at turn-sites, (c) Trajectory of the tangent 
vector on the unit sphere. An area of A(L) = 2tt is enclosed. 



of turn-sites keeping it within the plane (Fig.l). Then one defines A(L) to be 
the oriented area encircled in the right of the trajectory v(x). On the lattice 
Z 3 , A(L) takes values 

A(L) = ^tott (to e Z). (4) 

In two dimensions, the quantity A(L) takes values ±2ir and this signature corre- 
sponds to that of s ±1 in (2). Therefore (3) incorporates an essential ingredient 
of three-dimensional loops and is regarded as a natural generalization of (2). 

As is evident from the above construction, there is certain ambiguity for the 
value of A(L). First, because the trajectory is drawn on a closed surface of area 
Air, A(L) is well-defined up to An. Second, the signature of A(L) is changed 
when the picked direction of L is reversed. I require that n(-) in (3) absorbs 
this ambiguity. Hence, it should satisfy 

n(A) = n{A + An), (5) 
n(-A)=n(A). (6) 

Eqs. (4), (5), (6) imply that the fugacity function n(-) can be specified by five 
parameters n(A),A — 0, \n, n, |7T, 2tt. 

In spite of the above restriction on n(A), the model (3) includes many in- 
teresting cases. Consider, for example, fugacity 

n(A)^n S ( o L) (A), (7) 

with n G R and 



1 otherwise. 



The sum in (3) is then restricted to configurations which consist only of loops 
with the oriented area A = mod An. It should be interesting to compare the 
site entropy with that of the model with n(A) = n - It is also tempting to 
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Figure 2: Vertex weights TV for Z verte x[s](x 1 ). 



ask whether such an additional constraint changes the critical behavior or not. 
The present case reminds one of the fully packed loop model in two dimensions. 
Its universality class differs from that of densely packed loop phase when the 
additional constraint that the loop length must be even is imposed [12, 16]. 



3 Transfer matrices from local vertex weights 

In order to construct a local layer-to-layer transfer matrix for the loop model 
(3), I define a vertex model and show that it is equivalent with (3). 

The local degree of freedom z of the vertex model lives on each link (r, rie^), 
r € Z 3 . It takes one of three values <— , — », and —(empty). On each site, six 
neighboring link variables interact by the vertex weight W defined in Fig. 2, 
where s(lo) is a function that satisfies 

s{A 1 )xs{A 2 ) = s(A 1 + A 2 ) (9) 

and is specified further below. The partition function of the vertex model is 

Vertex [s] (x- 1 ) = II W({z((r,r±ei))}). (10) 

z=<-,->,- r£Z 3 

Evidently, the partition function (10) has a local transfer matrix which is a 
product of W's. 

Now I show that (10) for an appropriate s is equivalent with (3). Because the 
weight W is nonzero only when there is one incoming and one outgoing arrow, 
contribution to the partition sum (10) comes only from the set C of directed 
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loop configurations 



cec Lec(c) L reLnz 3 



n s ( w «) 



(ii) 



where s(w(r)) is the weight for the vertex at r e Z 3 in Fig. 2. 

The factor in the square bracket in (11) is associated with a directed loop 
component L. To evaluate this quantity, it is crucial to observe that 

A(L) = J2 "(r). (12) 

reLnz 3 

Actually, the weight system in Fig. 2 is designed to have this property in refs.[17, 
18, 19] in the context of random walk with a spin factor [20, 21]. Combining 
the properties (12) and (9), one finds that the factor in the square bracket is 
simply s(A(L)). 

Using the same trick as that used in (2), one can further write the partition 
function (11) as a sum over undirected loop configurations 

Z vcrtcx {s}(x- 1 ) = Y,x mc) - N II (s(A(L)) + s(-A(L))). (13) 

cGC LeC(c) 

Therefore, if n(A) in the loop model (3) can be written as 

n(A) = s(A) + s(-A) (14) 

with a function s that satisfies (9), then the vertex model partition function 
■^vertex [s] (x _1 ) defined above is equal to (3). 

The requirement (9) together with the restriction (4), (6) and (5) on n(A) 
forces s(A) to have a simple form 

s(A) = e lJA (15) 

with J e Z/2. It is enough [17] to consider the cases J = 0, |, 1, | and 2 because 
of (4). Hereafter, one introduces a shorthand notation 

Zjix- 1 ) = Z loop [n(A) = e UA + e"^]^" 1 ). (16) 

The vertex weights at J = and 2 enjoy a special property s(u>) = s(—u>). 
This enables one to define a model with only two microscopic states <-» and — : 

^ ertex W(x- 1 )= J2 II ^((r.riei))}) 

z=<->,- rGZ 3 

cec Lec(c) 
which I denote by Zqi and Zy . 
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Table 1 : The generalized fugacity n which generates the semigroup of allowed 
models. The function 5^ a \-) is defined in eq.(8). 

The fugacity functions corresponding toJ = 0,^,l,§,2,0' and 2' are listed 
in Table 1. They indeed give n(A) < or n(A) £ Q for some loops. 

Although only finite number of vertex models J = 0, |,2, 0' and 2' 
have been constructed above, it is possible to construct an infinite number of 
ones by taking the direct sum of the space of their microscopic states. More 
precisely, one generalizes the link variable to take one of 2q+ 1 (q E Z + ) states: 
^k,^k with the k labeling colors k = 1, . . . , q and an uncolored empty state 
— . Introducing parameters Jk € {0, |,1,|,2}, the vertex weight assignments 
in Fig. 2 are supplemented by additional rules: 

• If the both two arrows have the fc-th color, then W — e lJkUJ . 

• If the two colors do not agree, W = 0. 

The cases Jk = 0', 2' arc handled in the obvious manner. 

The fugacity of the 'direct sum' model Zj 1 $ ) j 2(S ... (S j q is simply the sum: 

n(A) = J2 [^ UkA + e ~" kA ) x B ^)\ > ( 18 ) 
fe=i 

1 ( J = 0, 2, 1, §, 2), 

I (j = o',n (iyj 

One immediately notices that 

Zo'®o' = %o, Z2'(£2' = Z2. (20) 

Thus the fugacity functions expressible via vertex models form an infinite semi- 
group under addition 2 generated by J = 0', 5, 1 5 § and 2'. 

2 This direct sum operation may be used for the lattice construction [17, 18, 19] of higher- 
spin three-dimensional field theories. 




G 



One can also take the 'tensor product' of the space of microscopic states of 
Zj 1 and Zj 2 to define a model Zj l8) j 2 . Let the link variable take five values 
(zi, z 2 ) =Tt, T4j IT; 11, an d II- The vertex weight W is defined to be the product 
of W's with J = Ji and J 2 . Then the loop fugacity becomes 

n{A) = [(e lJlA + e- lJlA ) x B{.h)} x [(e lJ * A + e - l - hA ) x B(J 2 )] . (21) 

However, a new fugacity function cannot be realized because the partition func- 
tion Z Jl ^j 2 is equivalent with an appropriate direct sum 

Zj 1 ®.j 2 = Z TH®-®l q ( 22 ) 
corresponding to the decomposition rule of the representation of SU(2). 



4 Entropy estimates 

I numerically diagonalize the transfer matrices constructed in section 3. Through- 
out this section, I concentrate on the fully packed limit x~ % = where all the 
sites are visited by a loop. This simple case is in fact very interesting case; in 
two dimensions, this limit yields a new universality class with a shifted cen- 
tral charge on several bipartite lattices and has been attracting much attention 
[16, 22, 23, 24]. It would be interesting to look at the limit where the two 
strong constraints are combined: the fully-packing constraint x^ 1 = and the 
constraint (7) on the shape of loops. 

The site entropy in the thermodynamic limit is defined by 

/[n](oo) = lim 1 \ogZ loop [n]{x~ l = 0)- (23) 

N—>oo 1\ 

I evaluate this quantity on quasi-one-dimcnsional geometry L± x L% x L%, — > 
oo while Li,L 2 kept finite by calculating the largest eigenvalue of the transfer 
matrix in an appropriate sector. 

Let T be the layer-to-layer transfer matrix in +e3 (vertical) direction for 
the vertex model defined in (10). Then T acts on linear combinations of arrays 
of Li x L 2 vertical (colored) arrows. One can take either hard-wall or periodic 
boundary condition in the horizontal directions. 

It is important to note that the transfer matrix T commutes with the oper- 
ator giving the net flow of arrows of fc-th color in +e 3 direction 

d k = (# Tfe) - (# ik), (24) 
which is understood as 

d k = (# Ifc) mod 2 (25) 
for Jfc = 0' and 2'. Thus T is block-diagonalized as 

T = 0T d , d =(<*!,... ,d q ). (26) 
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Table 2: The site entropy estimated numerically. L± x L 2 is the size of a layer 
while (p) and (h) mean periodic and hard-wall boundary conditions in a layer. 



The quantity (23) is obtained as 

/[n](oo)= lim J—l og \\° (L u L 2 )\, (27) 

where \ l A {L\, L 2 ) is the «-th largest eigenvalue of T^{L\,L 2 ). The condition d = 
excludes unwanted configurations that have unbalanced arrows traveling along 
the infinite direction. Shown in Tables 2 and 3 are the finite- Li, L 2 results 3 . The 
asymmetric Lanczos algorithm is utilized for the present sparse eigenproblcm. 

The obstacle in making L\ and L 2 large in the actual numerical work is 
of course the exponential growth of the dimensionality of the transfer matrix. 
The selection of d = sector helps to reduce the dimensionality, though the 
improvement is polynomial. For a direct sum Zj l(S ... (S j q , the dimensionality is 
estimated to be 

/ 3 \ "/ 2 

dim(d = sector) ~ (1 + 2u + p) LlL2 x 2~ p - - - , (28) 

\4ttLiL 2 ) 

where p = #{fc| J k = 0' or J k = 2'}, u = q — p. One notices that it is better to 
recast Z 2 i^ 2 i as Z 2 if one is interested in its d = sector. 

The exponential growth is severe even after restricting to the d = sector. 
In order to increase L\L 2 as much as possible within the available computer 
resources, I have further decomposed T with respect to the eigenvalue of shift 
(lattice momentum) operator for the periodic boundary case where the transla- 
tional symmetry is present. I have looked at the zero-momentum sector Jj/ ' ^ 
as described in the appendix A. It is quite natural to expect that the largest 

3 I have also measured several leading eigenvalues of with d = (di, . . . , d q ), = 0, ±1. 
These are related to correlation length of operators in the theory. These results will be reported 
elsewhere. 
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Table 3: The site entropy estimated numerically. L\ x L-i is the size of a layer 
while (p) and (h) mean periodic and hard-wall boundary conditions in a layer. 

eigenvalue lies there. By this decomposition, the dimensionality of the eigen- 
problem is reduced, at most, by {L\L^)~ X . As a drawback, the matrix Ji/ 0,0 ) 
becomes less sparse than the original To- With both effects combined, some 
improvements in the memory usage and the CPU time are observed. Thus the 
analysis of larger systems becomes possible for the periodic boundary 
seen in Tables 2 and 3. 

5 Discussions 

In the comparison between the periodic and the hard- wall boundary conditions, 
one notices that the periodic case has always larger site entropy. This is because 
many of the loops that wind non-trivially in the horizontal directions satisfy 
A{L) = and the loops with A(L) = contribute to every partition sum with 
a positive fugacity. 

The numerical works in the present study have been carried out on modest 
workstations. Unfortunately, information in the thermodynamic limit L±, L2 —> 
00 is out of reach in the present analysis. For the study of criticality, ref.[17], 
where random walks with the weight in Fig. 2 are studied, is quite suggestive. 
It is reported that Euclidean symmetry is not always recovered even in the 
continuum limit. 

I discuss relations with combinatorial problems below. 
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5.1 Even and odd number of loops of a specific type 

For most allowed values of J, the loop fugacity takes both positive and negative 
values. Some interesting combinatorial information is encoded in these mod- 
els. For example, the linear combinations \{Z i ± Z 2 >) counts the number of 
loop configurations such that there are even (odd) number of loops for which 
§A(L) = 1 mod 2, e.g., 

-{Zv±Z 2l ) = { — * (29) 

where C ovon ( odd ). a is the subset of C and is defined by the following properties. 

• c e C even ( dd),o contains even (odd) number of loops with A{L)/aiT = 1 
mod 2. 

• All other loops in c e C von(odd).a satisfy A(L)/an = mod 2. 

Similarly, the quantities ^(i?o'©2' and ^(Zo'§ 2 '§i ± -Zi^a) are interpreted 

as the sums over C ovon (odd),i and C ovon ( dd),2- 

Of the two Z's in (29), the one with the larger leading eigenvalue dominates 
the sum in the limit L 3 —> oo studied in section 3. In finite geometries, both 
terms contribute to yield an exact number. 

5.2 Self-avoiding walk 

The partition function (1) in the limit n — > corresponds to the enumeration of 
self-avoiding walks. Self-avoiding walks in three-dimensions have mainly been 
studied by the exact enumeration method due to the lack of transfer matrix 
formalism as pointed out in the introduction. 

The model I propose in this paper can be regarded as a step forward to 
overcome this difficulty; in the models i?o'e2'> i?o'©2'©i an d ■^o'©2'e i ei©-' ^ ne 
fugacity is set zero for families of loops. This is, however, achieved by paying 
the cost of having larger n for another family of loops. Within the present 
construction, loops with A(L) = mod 4ir cannot have weights different from 
the number of possible link states. Thus the partition function listed above 
serve only as a very loose upper bound for the entropy of self-avoiding walks. 

The problem of construction of a local transfer matrix to enumerate self- 
avoiding walks on three-dimensional lattices still remains open. 

5.3 Mapping to ribbon configurations 

The oriented area defined in (3) has a nice geometric interpretation as holon- 
omy. The tangent vector v(x) moves along a trajectory on S 2 . Let the unit 
vector tangential to this trajectory at v(x ) be u (Fig. 3). Consider the paral- 
lel transport (in the sense of riemannian geometry) of u along the trajectory 
v(x) on S 2 . 
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Figure 3: When uo, the tangent vector at v(xo), is parallel-transported along 
the trajectory on S 2 , it receives holonomy whose magnitude (angle 9) is equal 
to the oriented area A(L) the path encloses. 



When u is transported back to v(x ), it gains some holonomy ( the angle 9 
in Fig. 3). This holonomy angle is given by the integration of the scalar curvature 
of S 2 over the domain encircled by the trajectory and is nothing but the oriented 
area A(L). In the real space Z 3 C R 3 , the holonomy described above is nicely 
kept track of by broadening the loop segment to a 'ribbon' with the distinction 
of the right and the reversed sides. The parallel transportation can be recasted 
as a rule of bending ribbons on sites, which is shown in Fig. 4. Among the 
loop configurations shown in Fig. 4, the partition sum Z\ oop with fugacity n = 

Sq , 5q ' and S 2 receives contribution from {(a)}, {(a),(b)}, and {(a),(b),(c)}, 
respectively. 

(2) 

For n — 4<5q , the sum in (3) is over ribbon loop configurations without mis- 
match. In that interpretation, the coefficient four is naturally regardes as the 
number of directions the right side can face. Therefore, Z\ oop [n = 4(5q 2 '] (a; -1 ) is 
nothing but the generating function of the number of allowed ribbon configura- 
tions 4 . 

Similarly, the partition sum for n = 25^ can be interpreted as the sum over 
the configurations of ribbons without the distinction between the right and the 

(I) 

reverse sides, while in the case n = 5 2 , the loop segment is just a chord. This 
interpretation suggests that (3) may be regarded as a model of polymers with 
various partially broken axial symmetry by, for instance, the presence of side 
chains. 

4 This observation clearly indicate a simple way of constructing transfer matrices for this 
fugacity function. The link variable represents the ribbon which can be placed in four ways. 
The rule of bending is implemented in the vertex weight. The construction in section 3, 
however, has an advantage; the size of the matrix can be reduced more by the choice of the 
sector d = 0. 
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(a) (b) (c) 

Figure 4: Examples of allowed ways of bending ribbons at sites (except for the 
ones indicated by gray arrows) The holonomy is accumulated at the sites indi- 
cated by gray arrows. The right and the reverse sides of ribbons are represented 
by white and black colors. 



5.4 Comparison with the connectivity basis 

The connectivity basis[25, 26] is very powerful in that one can always write a 
transfer matrix for a loop model with respect to it. It has been very useful for 
numerical calculation in two-dimensions. 

Nevertheless, I have avoided the use of the connectivity basis in the present 
work. The reason is the following. First, its fundamental degrees of freedom 
are not the link variables and the transfer matrix with respect to it is not local. 
Local transfer matrices have merits even in two dimensions. Namely, it paved 
the way to the Bethe ansatz solution[27, 28] and the conformal field theoretic 
description [29, 30] via Coulomb gas representation. Second, in three dimensions 
and higher, the size of connectivity basis grows considerably because of the lack 
of the planarity constraint. It is not clear if it is effective to perform numerical 
calculation in this basis. In two dimensions, the present basis is as good as the 
connectivity one [12, 22]. 

I suppose it is very important to see how useful the connectivity basis in 
three dimensions is and to try to improve the efficiency of the calculation in 
that basis. 
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A Projection to zero momentum subspace 

In this short note, I describe the block-diagonalization of the transfer matrix 
with respect to the eigenvalues of the lattice momentum operators. The zero- 
momentum subspace and a reduced transfer matrix which acts on it are explic- 
itly constructed 5 . 

One may start with the (2q + l) LlL2 -dimensional whole space of colored 
arrow configurations or an eigenspace of the operator d. One considers the 
matrix elements in the basis Uj, (i = 1, ... ,m), each of which represents a 
single arrow configuration such as jj. | • • • j | J.: 

m 

In this natural basis, the matrix T becomes sparse. 

Let S\ and S2 be discrete shift operators in the horizontal directions. Then 
the vectors 

Li-l L 2 -l 

v i = E E (Si) a (S2)\ 

a=0 6=0 

are zero-momentum ones. 

One classifies the index set as {1, . . . , to} = Ujf^Vj by an equivalence rela- 
tion i ~ j Vi = Vj. Then 1 = 1,... ,M labels the zero- momentum subspace. 
The (I, J)-component of the block matrix is simply 

( To (°.°>)/ = T i j (* e Vj) 

This procedure is fairy easy to implement in the sparse algorithm. 

Evidently, a slight modification of the above procedure enables one to focus 
on a chosen non-zero momentum subspace. It will be useful for identifying 
excited states. 

It is noted that the above block-decomposition can be applied even if the 
seam factor is present, e.g. to two-dimensional O(n) model with cylinder topol- 
ogy. One can make the system translationally invariant by distributing the seam 
factor among all horizontal links. I have checked that this prescription improves 
the efficiency of the the enumeration of Hamiltonian cycles perfomed in ref. [22] 
although the weight system becomes system size dependent. 
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